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ABSTRACT 



Flow rate measuring method and apparatus for deter- 
mining a flow rate of a fluid flowing through a pipe by 
using a fuzzy inference system which is operative on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 
and membership functions defining meanings of state- 
ments described in said IF and THEN parts. Row ve- 
locity of the fluid is measured at a plurality of locations 
within the pipe by a plurality of sensors. A table of 
membership functions indicating grades of fitness of the 
flow velocities to a distorted flow model distribution is 
provided in a storage together with a mean flow rate 
calculating expression for deriving a flow rate corre- 
sponding to the distorted flow model distribution. A 
membership grade of the measured flow velocity to the 
distorted flow model distribution is determined by re- 
ferring to the membership function table through the 
medium of the measured flow velocity, and a weighted 
mean of the flow rates given by the mean flow rate 
calculating expression is determined by using as a 
weight the membership grade to the distorted flow 
model distribution. 
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METHOD AND SYSTEM FOR MEASURING 
FLUID FLOW RATE BY USING FUZZY 
INFERENCE 

BACKGROUND OF THE INVENTION 

The present invention relates to a method and a sys- 
tem for measuring a flow rate of a fluid flowing within 
and through a pipe. More particularly, the invention is 
concerned with a flow rate measuring method and a 
system therefor with are capable of determining or 
measuring stably a mean flow rate of a fluid flowing 
through a pipe such as tube, conduit or the like in which 
distorted flows and disturbances such as turbulences 
may exist. 

For measurement of flow rate of a fluid flowing 
within a pipe, there has heretofore been known a plural- 
sensor type flow meter system in which a plurality of 
detection elements (hereinafter referred to as the sen- 
sor) are disposed within the pipe and a single-sensor 
type flow meter system in which' only one sensor is 
installed within the pipe. 

As a typical plural-sensor type flow meter, there may 
be mentioned those disclosed in JU-A-6 1-1954188 and 
JP-A-6M20016, respectively. The first mentioned flow 
meter is based on the concept that the output signals of 
the individual sensors are integrated for determining 
stably an approximately true means (average) flow rate 
by using a calculating circuit of a corresponding struc- 
ture notwithstanding of presence of distortions or devi- 
ations in the flow velocity (i.e. a distorted flow velocity 
distribution) brought about by geometrical factors (e.g. 
shapes) of an upstream pipe section, as is disclosed in 
JU-A-61-195418. 

In the single-sensor type flow meter, a single hot wire 
element is installed within the pipe under tension in 
order to determine a curvilinear integral (line integral) 
of a flow velocity for thereby obtain a mean (average) 
flow rate, as disclosed in JP-A-54- 125060. 

The plural-sensor type flow meter in which a plural- 
ity of sensors are installed within the pipe suffers from a 
problem that when any one of the sensors is located in 
a region where flow disturbances due to vortexes, flu- 
idal separations or the like brought about by geometry 
or shape of the pipe are significant, infinitesimal varia- 
tions or fluctuations (fluidal noise) of the flow velocity 
aTe directly fetched by the mean flow calculating circuit 
as the input thereto, making it difficult or impossible to 
obtain the mean flow rate in a stable manner, because 
weights imparted to the outputs of the individual sen- 
sors are selected to be equal to one another or Axed. 

The single-sensor type flow meter mentioned above 
also suffers the similar problem because no measures are 
taken against the influence of the infinitesimal variations 
in the flow velocity. 

It is further noted that the location where the infini- 
tesimal variation takes place in the fluid flow velocity is 
movable as a function of change in the distorted flow 
distribution profile or pattern. Consequently, the mean 
flow rate calculation using the fixed weight can not 
stably or consistently assure the availableness of accu- 
rate flow rate value by flexibly coping with the infinites- 
imal variation, giving rise to another problem. 

SUMMARY OF THE INVENTION 

It is an object of the present invention to provide a 
flow rate measuring method and a system for carrying 
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out the same which can evade the shortcomings and 
problems Of the techniques mentioned above. 

Another object of the present invention is to provide 
a method and a system for measuring a flow rate which 
are capable of measuring a true mean flow rate even in 
the presence of distorted flows and disturbances by 
using a fuzzy inference and a fuzzy integration in deter- 
mination and calculation of a mean flow rate of a fluid 
flowing through a pipe on the basis of outputs of plural 
sensors installed within the pipe. 

In view of the above and other objects which will be 
more apparent as the description proceeds, there is 
provided according to a first aspect of the present in- 
vention a flow rate measuring method of determining a 
flow rate of a fluid flowing through a pipe by using a 
fuzzy inference system which comprises the steps of 
measuring flow velocities of the fluid at a plurality of 
locations within the pipe, preparing a table of member- 
ship functions indicating grades of fitness of the flow 
velocities to a distorted flow model distribution, prepar- 
ing a mean flow rate calculating expression for deriving 
a flow rate corresponding to the distorted flow model 
distribution, determining a membership grade of the 
measured flow velocity to the distorted flow model 
distribution by referring to the membership function 
table through the medium of the measured flow veloc- 
ity, and calculating a weighted mean of the flow rates 
given by the mean flow rate calculating expression by 
using as a weight the membership grade to the distorted 
flow model distribution. 

According to another aspect of the invention, there is 
provided a flow rate measuring system for determining 
a flow rate of a fluid flowing through a pipe by using a 
fuzzy inference system which comprises a plurality of 
flow sensors for measuring flow velocities of the fluid at 
a plurality of locations within the pipe, a storage con- 
taining a first membership function table indicating 
grades of fitness of a mean value of the flow velocities 
within the pipe to a plurality of distorted flow velocity 
levels, respectively, a storage containing a second mem- 
bership function table indicating grades of fitness of the 
flow velocities to a plurality of distorted flow model 
distributions for the flow velocity levels, respectively, a 
storage containing mean flow rate calculating expres- 
sions for deriving flow rates corresponding to the dis- 
torted flow model distributions at the flow velocity 
levels, respectively, a unit for determining a mean value 
of the flow velocities measured at the plurality of loca- 
tions within the pipe to thereby determine on the basis 
of the mean value first membership grades of the mean 
value to the flow levels, respectively, by reference to 
the first membership function table through the medium 
of the mean value, a unit for determining second mem- 
bership grades of the flow velocity values to the dis- 
torted flow model distributions for the flow velocity 
levels, respectively, by reference to the second member- 
ship function table through the medium of the measured 
flow velocity values, and a unit for calculating a 
weighted mean of the flow rates given by the mean flow 
rate calculating expressions, respectively, by using as 
weights the first and second membership grades, for 
thereby determining the flow rate of the fluid within the 
pipe. 

The present invention starts from the facts that the 
flow velocities at various points in a fluid flowing 
through a pipe exhibit a frequency distribution having a 
unique mean (average) value and variance due to influ- 
ences of distorted flows and disturbances and that the 
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frequency distribution profile is characteristically inner- FIGS. 6A to 6C are views showing frequency distri- 

ent to the geometries of the pipe line and the mean flow butions at different flow velocity levels for illustrating 

velocity prevailing within the pipe. the contents of the table shown in FIG. 5; 

For calculating the mean flow rate by taking advan- FIG* 7 is a view showing a structure of a distorted 

tage of the characteristics mentioned above, it is pro- 5 flow model membership function table; 

posed as an embodiment of the invention that the de- FIGS. 8A to 8C show three exemplary configura- 

tected values outputted from a plurality of sensors dis- ^ ons of sections through which a fluid under mea- 

posed within the pipe are made use of as fuzzy variables, surement flows; 

wherein software and hardware are so structurized that ? GS - 10 * J « Y 1 ^ £ 0Wm S frequency distri- 
correspondences are established by the fuzzy inference 10 but,ons °J fiow , velocities in distorted flow model pat- 
between or among the fuzzy variables, a membership te ™J? dependence on particular shapes of pipes; 
function table (first membership grade function table) , JF 10 - 10 shows a «™^ c of a mcan flow ratc adcu ' 
containing information about the levels to which the latmg expression toble l^; 

flow velocities belong, a distorted flow model function 1€ FI °* 11 15 8 fl ™ chart for grating an inference 

table (second membership grade function table) con- 15 J™*™ executed by a mean flow rate mference part 

taining distorted flow distribution model patterns and nG .n'shows a structure of the flow velocity level 

an empirical expression (mean flow rate calculating mcmbcrshi funclion toUe for illustrating operation of 

expression) for determining mean flow rates corre- the flow ^ measuring system according t0 \ second 

spending to the flow velocity levels and the distorted 20 embodimcnt of ^ invention; 

flo w distribution patterns, respecuvely FIG. 13 is a view showing a structure of the memory 

When a distorted flow distribution of a given pattern emp i oycd ^ thc flo w rate measuring system according 

exists in the vicinity of plural sensors, the fitness grades t0 a third embodiment of the invention; 

(similarity grades) or membership grades of the dis- 14 is a flow chart for illustrating the inference 

torted flow distribution to a plurality of distorted flow 2 $ procedure executed by a mean flow rate fuzzy inference 

distribution patterns prepared previously independent part 0 f tne system according to the second embodiment 

of the former are determined by fuzzy arithmetic opera- of the invention; 

tion, whereon the mean flow rates which correspond to FIG. 15 is a view showing another example of instal- 

the distorted flow model patterns, respectively, are lation of sensors within a pipe; 

determined in accordance with the empirical expres- 30 FIG. 16 is a view showing a structure of a fuzzy 

sions corresponding to the distorted flow model pat- arithmetic unit according to a fourth embodiment of the 

terns. Further, by using the membership grades (fitness invention; and 

grades) as weight, a weighted mean value is calculated FIG. 17 is a view showing a general arrangement of 

for the distorted flow model patterns. an air suction system in an internal combustion engine 

By determining the membership grades for the dis- 35 to which an embodiment of the flow rate measuring 

torted flow model patterns, as mentioned above, it is system according to the invention is applied, 

also possible to calculate mean (average) flow rates for DESCRIPTION OF THE PREFERRED 

those distorted flow distributions which are deviated EMBODIMENTS 
from the previously prepared distorted flow model 

patterns. In this manner, the capability of the flow rate «0 Before entering into detailed description of preferred 

measuring system to cope with variations of the dis- embodiments of the flow rate measuring method and 

torted flow distribution can significantly be enhanced. svstem according to the present invention, elucidation 

Further, since the mean flow rate is not straightfor- will first be made of the principle or concept underlying 

wardly determined by using directly the output values a l yP Ical example of the mean flow rate calculation 

of the sensors, the mean flow rate as calculated is essen- 45 i* u 8 ht bv _ invention, 

tially insusceptible to direct influence of local disturb- FIG r 2 * 18 8 diagram for illustrating gener- 

ances taking place around particular sensors. Thus, ation of distorted flows in a fluid flowing through a 

there can be realized a flow rate measuring system lri . 2A, a reference character A denotes gen- 

which is highly stabilized against disturbances of the <n "f* * A ? 1V *\u"* ^ 

fluid flow under measurement. 50 ^f"* t0 T L^™ 8 * 1 ^ L^^?^ 

under measurement flows, and B denotes distorted 

BRIEF DESCRIPTION OF THE DRAWINGS flows of the fluid brought about within the pipe A under 

, - . uu i, j- v v ii the influence of geometrical factors such as shape of a 

FIG 1 is a block dtagram showing schematically a . ^ Furth a fef 7 rence 

general arrangement of a flow rate measuring system „ ch VracterC denotes generally a flowing direction of the 
according to a first embodiment of the present uiven- fluid Now> lcf$ ^ a flow velocity measurement 

t,0 JH^ • . + .„ ... of a fluid flowing through the pipe line A. So long as the 

FIG 2A is a schematic diagram for illustrating dis- ^ flow rate of the fluid ^ pjpc Une A is 

t0 ~l fl ™ S ° f 3 . nuid . floWUlg thrOUgh a pipe; constant, the ideal conditions for the measurement of 

FIG. 2B is a view showing graphically a frequency « flow velocities at measuring points Pi, P 2 and P 3 within 

distribution of distorted flows; the pipc A should be such that flow velocity values ui, 

FIG. 3 is a flow chart for illustrating a typical exam- U2 and u 3 measured at the respective measuring points 

pie of procedure for calculating a mean flow rate ac- p 1( p 2 and P 3 remain always constant with variance 

cording to the teaching of the present invention; among the frequency distributions of the flow velocity 

FIG. 4 is a view showing an internal structure of a 65 measurement values ui, u 2 and m being zero. However, 

memory used for fuzzy inference of the mean flow rate; in practical flow velocity measurements, the fluid flow 

FIG. 5 is a view showing a structure of a flow veloc- is inevitably accompanied with deviations and infinites- 

ity level membership function table; imal deviations under the influence of disturbances such 
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as vortexes, fluidal separation and/or the like phenora- Referring to FIG. 3, execution of the procedure starts 
ena taking place upstream of the measuring points as from a step SO. In a step SI, detected or measured flow 
well as under the influence of the geometrical configu- velocity values ui, U2, . . . , uw are fetched from flow 
ration of a pipe section located upstream. As a conse- sensors disposed at the measuring points Pjm . . . , 
quence, the actual flow velocity (vj, V2, V3) exhibits 5 p^, respectively. Subsequently^ in a step S2, a flow 
such frequency distributions which differ from_ one velocity mean (average) value u of ui, U2, . . . » mm is 
another in respect to the temporal mean values \\ f V2 determined on the basis of the flow velocity measure- 

and V3 and variance, as illustrated in FIG. 2B. In this me nt values ui, U2 mm in accordance with 

conjunction, it is noted that differences among the tem- 
poral mean values of the flow velocities vi, \i and V3 10 ^ m 
characteristically represent distribution of distorted u = 1 tt/M 
flows, while variance at the individual measuring points . . . . 
character^ the deviation or fluctuation of the flow Next ' m a ste P S ?/ a flow measurement value 
velocity brought about by the disturbances or turbu- range corresponding to a flow veloaty classification 
lences such as vortexes, fluidal separation and/or the 15 ™?gf ^ which the mean velocity of the fluid flow 
like phenomena. When the mean flow rate within the w,tlun the P l I* A bclon « s and a ^stoned flow model 
pipe is to be determined on the basis of a plurality of the membership function table bearing a strongest or best 
flow velocity measurement values u h u 2 and u 3 exhibit- correspondence relation to the mean flow velocity are 
ing the distributions described above, it may be con- determined on the basis of the abovementioned mean 
ceived to adopt a method of determining a simple mean 20 fiow velocity value u in accordance with the following 
value of the measured flow velocities u\ f u: and U3 or a expression (2). 

method of determining a weighted flow rate mean value J«Klfo)7KEWM« 0) 

(wivi+w2V2-|-w3V3) by weighting the measured flow 
velocities ui, u 2 and 113 with coefficients wi, W2 and W3, where 

respectively, which can be derived from Reynolds num- 25 ^presents a constant m a range of 1, . . . , K, 
ber (Re) representing an approximate value of the mean U " *F™?** a fl ° w ve,0 f ltv classification range to 
flow rate by taking into account the geometrical config- t wh,ch the mean fiow velocity vdue u belongs and 
uration of the flow path and other factors. However, 1 ^presents a specified or particular range of flow 
the simple mean method is disadvantageous in that flue- velocity measurement values, 

tuations or small excursions will inevitably be included 30 lli order t0 ™ ow ^ distributions of distorted flows 
admixefiy in the mean value as determined, because no 10 be handled " fuzz V variables, the flow velocity mea- 

consideration is paid to variance of the flow velocity surement values ui, u 2 UA#are divided by a maxi- 

due to the infinitesimal variation of fluctuation of the mum flow velocity value u max given by 
flow velocity. On the other hand, the weighted mean 

method will encounter difficulty in determining defi- 35 «*.... «w) (3) 

nitely the distribution of distorted flows because of the A A « , . .. , - . - 

variance of the How velocity, which in tun, means that t0 there * y < * t "» "onnahzed flow velocity nieasure- 
it becomes difficult to determine evaluation coefficients m "V ValueS <?«• V 2 U *] ^ j* md S % t . 
related to the distribution profile or pattern of the dis- . Subsequently, m a step S6, i. «jletermmed on the 
torted flows which are to be taken into consideration in 40 "ft 0 "? ^""n f* bdow 

determinating the applicable or pertinent weights w,, ° f the dlst . orted flow modd J" 11 "™, ±e "^ized 
W2and W3 flow velocity measurement values uj f U2, . . . , upbears 

At this juncture, there may be mentioned two main ? e best simUarity of fitness (membership grade) The N 
factors which make it difficult to calculate the mean ^zy inference rules are prepared m correspondence to 
flow velocity They are* 45 SP* 0 ™ flow velocity levels. More specifically, N dis- 

(1) fuzziness of the mean flow rate brought about by ^ floW , IDOdcl patternS m combined with plural 
variations in the mean flow rate level within the pipe, fl ° w r veloc I tv measurement values in the form of a 
gjjjj table (refer to, for example, FIG. 7), whereon member- 

(2) fuzziness of the mean flow rate due to variance in or fltness « rades (membership function val- 
the flow velocity measurement values due to the dis- 50 ues) are allocated to the abovementioned combinations, 
torted flow distribution and the disturbances such as respectively. 

mentioned above. Expressing mathematically, 

Accordingly, by resorting to the use of a fuzzy infer- 
ence method for determining the mean flow velocity ^ tE Q t VjU 
level and the profile or pattern of the distorted flow 55. 
distribution to thereby take into account the fuzziness in 92 Vnk * 
determination of the weighting coefficients mentioned IF 
above, it is believed to make steady (or stabilize) the 
. calculation of the average or mean flow rate, as will be *m -EQ- *W*) 
described below in detail. 60 

Representing in general by M the number of points at THEN 
which the fluid flow velocity is to be measured, descrip- *=fyfaif?h »a/) <*) 

tion will be made on a method of determining the mean 

flow velocity level by resorting to a fuzzy inference where i=l, 2, , . . , N, and Jfe=l,2 K. 

method. FIG. 3 is a flow chart illustrating a typical 65 In the above expression (4), vi, \i> ...» vm are fuzzy 
procedure for calculating a mean fluid flow rate by variables representing typically the flow velocity mea- 
resorting to a fuzzy inference method according to an surement values obtained at the sensor locations or 
embodiment of the invention. measuring points Pi, P2, . . . , Pjv within the pipe, which 
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are normalized by the method described previously, and rate calculating function fa are averaged by weighting 
V/ut, Vfljt, . . . t Vim are fuzzy variables representing the values with the similarity or fitness grade (member- 
the fuzziness of the flow velocity distributions at the ship grade) 4>s* for the i-th distorted flow model pattern 
abovementioned sensor locations Pi, P2, . . . . Pjif, re- (p mi ) of the flow velocity measurement values (ui, U2, . 
spectively, which are derived from the frequency distri- 5 . . , uj*) at the k-th flow velocity level U*. 
butions of flow velocities at the sensor locations at a Next, on the basis of the mean flow rate uo* deter- 
given level of the mean flow velocity within the pipe mined as described above, a mean value uo is determined 
and normalized through division by a maximum distri- m accordance with the undermentioned expression (7) 
bution frequency. Further, f/jt represents a functional by using as the weight the membership grade {« for the 
expression for calculating the mean flow rate within the 10 k-th flow level determined previously (a step S8J, 
pipe on the basis of the flow velocity measurement whereon the calculation processing comes to an end 
values (ui, U2, . . . t UA/)» wherein the functional expres- ( stcp 
sion may include empirical rules or expressions. More 

specifically, f# represents an expression for calculating ^ ^ 

the mean flow rate corresponding to the maximum 15 1 ^mk 

distribution frequency in the distorted flow distributions « 0 = 

as the i-th fuzzy inference rule. 2 tki 

In the expression (4), "IP* part is referred to as the 

condition part with 'THEN" part as the conclusion M . _ . , 

part. Notation EQ." does noi represent equality of 20 . » germination of the mean ^w rate in |»«dnfle 
values but represent an arithmetic operation for check- with ^ above expression (7), the^ weights for he 
ing the membership grade of the fuzzy variable v to the individual flow velocity levels and the weights for the 
fuzzy variable V, which will be described below in membership grades to the mdividual distorted flow 
dcla il model patterns are synthetically taken into consider- 

It is now assumed that the flow velocity measurement 25 ation, which in turn means that inference of the mean 
value ui detected at the measuring point Pi is outputted flow rates for intermediate distorted flow distributions 
at a given time point In that case, the membership which are not contained in the distorted flow model 
function indicating to what extent the flow velocity pattern classification table can be realized as well. Fur- 
measurement value ui contributes to validation of the ther, by virtue of the nature of the fuzzy inferences 
condition part, i.e. the grade of equality or fitness of the 30 given by the expressions (4) and (5), unwanted influence 
value ui to the fuzzy variable V,* is represented by to the inferred value of the mean flow rate can be initi- 

al), where ui represents the normalized flow velocity gated even when any one of the flow velocity measure- 
measurement value at the measuring point Pi. Further, ment values exhibits significant fluctuations due to the 
the membership function of the condition part for the disturbances, because the membership grade (weight) 
flow velocity measurement value (ui, U2, . . . , ujy) is 35 <f>,is decreased correspondingly, 
defined as follows: Now, the present invention will be described in detail 

and concrete in conjunction with preferred or exem- 
(5) plary embodiments thereof. 
_ ^ FIG. 1 is a block diagram showing schematically a 

"yQ^ 40 general arrangement of a flow rate measuring system 

according to a first embodiment of the invention. In this 
where a symbol * W represents an operand indicating figure, a reference numeral 1 denotes a fuzzy arithmetic 
that the membership function assumes a minimum (operation) unit which incorporates therein a processor 
value. On these conditions, it is possible to determine unit 11 and a memory 12, wherein the processor unit 11 
with the expression (4) to what extent the flow velocity 45 and the memory 12 are interconnected through an inter- 
distribution of the measurement values (uj, U2, . . . , ua/) face 30 so as to allow data transfer therebetween. The 
bears "similarity or fitness" to the distorted flow model memory 12 may be constituted by a random access 
pattern classified as the i-th distorted flow model pat- memory or a read-only memory. On the other hand, the 
tern, when judged on the basis of M flow velocity mea- processor unit 11 includes a central processing unit 
surement values ui, . . . , um, wherein $uc given by the 50 (CPU) 110, a memory 112, an input/output (I/O) circuit 
expression (5) indicates the extent or grade of "similar- 114 and a bus 116 and is in charge of overall control of 
ity or fitness", i.e. the membership grade. In this con- the whole system. Further, the processor unit 11 is 
nection, a minimum value of the term \iijk appearing in provided with a data bus 1110 for outputting externally 
the expression (5) means that the similarity or member- the results of calculations. Additionally, there are pro- 
ship grade . is evaluated on the severer conditions. 55 vided for the processor unit 11 input terminals 1111, 
Through similar procedure, the membership grades 4>tk H12, . . . , 111M which are connected to plural (M in 
are determined for N inference rules, respectively, the illustrated embodiment) sensors 51, 52, ... , 5M, 
whereon the mean flow rate within the pipe is calcu- respectively, by way of an analogue-to-digital (A/D) 
lated in accordance with the following expression: converter 2, a noise eliminating filter 3 and respective 

60 detection circuits 41, 42, .... 4M so that the processor 
N (6) unit 11 can fetch data from the sensors 51, 52, ... , 5M, 
1 Matou n ■»*) respectively. The sensors 51, 52 5M are disposed 



A r within a pipe 100 through which a fluid flows at the 

it 1 locations appropriate for detecting distorted flows gen- 

65 eralty denoted by a reference numeral 101. In this con- 
where k= 1, 2, . . . , K. junction, it should be mentioned that the sensors 51, 52, 
The above expression (6) equivalently means that . . . , 5M may be disposed on a same plane as shown in 
values determined in accordance with the mean flow FIG. 1 or on mutually different planes as illustrated in 
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FIG. 15 on the condition that these sensors are posi- FIGS. 6A to 6C are views for illustrating, by way of 

tioned at different heights h. example, the contents of the flow velocity membership 

The sensors 51, 52, .... 5M are constituted by an function table shown in FIG. 5. In these figures, there 
element capable of detecting a change in a physical are graphically illustrated the flow velocity member- 
quantity (hereinafter referred to as the physical change) 5 ship functions in three different flow velocity ranges 
brought about by motion or kinematic behavior of the Vdu U<n and U& into which the mean values u of the 
fluid in which the sensors are immersed and may each flow velocity measurement values ui, . . . , u m derived 
be constituted by a pressure sensor, a thermistor, a hot from the outputs of the sensors 51, 52, ... , 5M are 
wire anemometer or the like. The detection circuits 41, classified (L thus being equal tojhree). Tbey are a low 
42, . . . , 4M serve for function to convert the physical 10 flow velocity range Vd\ Ou^u^Urfi), a medium flow 
changes detected by the sensors 51, 52, . . . , 5M into velocity range XJ<n <jy i <u ^ \idi) and a high flow veloc- 
electric signals, respectively. To this end, detection ity range Urf3 (u<e <u«Urfj). Referring to FIGS. 6A, 6B 
circuit may be constituted by a commercially available and 6C taken along the abscissae are three flow veloc- 
unit selected to conform with the characteristics of the ity levels Ui, U2 and U3, respectively, into which the 
sensor. It should further be mentioned that the sensors 15 mean flow velocity u is classified (with k-3), i.e. a low 
51, 52, . . . , 5M may be of such type as to operate on flow velocity level UiJuag^u^Uai), a medium flow 
mutually different principles with the detection circuits velocity level U2 (u«i <u^ iia2) and a high flow velocity 
41, 42, , . . , 4M being then adapted to operate on corre- level U3 (u tf 2<u«Ua3). while taken along the ordinates 
spondingly different principles. the membership values (fitness or membership grades) 

The signals detected by the detection circuits 41, 42, 20 {*/of the flow velocity levels Ui, U2 and U3 in the flow 

. . . , 4M are inputted to the noise, eliminating filter 3 velocity ranges Vdh V<n and U<«» respectively. It will 

through which high frequency noise components such be appreciated that the flow velocity level membership 

as electric noise, thermal noise and others are elimi- function table shown in FIG. 5 contains the member- 

nated. Output signals from the filter 3 are converted ship grades (fitness grades) {jt/of the flow velocity lev- 

into digital signals by the A/D converter 2 to be subse- 25 els shown in FIGS, 6A, 6B and 6C in the flow velocity 

quently supplied to the fuzzy arithmetic unit 1. ranges Vd\ to Vdt (L=3 in this case), respectively. In 

Through the signal transmission path mentioned this manner, when the flow velocity mean value u be- 
above, the flow signals detected by the sensors 51, 52, . longs to or falls within the low flow velocity range Vdu 
. . , 5M are applied to the input terminals 1111, 1112, . . the membership grades £1 u li\ and (31 of the flow ve- 
. , 111M as the signals which can discriminatively be 30 locity levels Uj, U2and U3 can be determined to be, for 
identified from one another. example, 1.0, 0.5 and 0.2, respectively, from the table 

Implemented in the memory 12 are a mean (average) content shown in FIG. 6A. Parenthetically, data illus- 

flow rate inference part incorporating inference proce- trated in FIGS. 6A to 6C are only for the illustrative 

dure and tables which are required for the processor 11 purpose. It will be understood that L and K may assume 

to perform fuzzy inference of the mean flow rate within 35 given positive integers, respectively, 

the pipe A on the basis of the signals applied to the input FIG. 7 shows a typical structure of the distorted flow 

terminals 1111, 1112, . . . , 111M. FIG. 4 is a view show- model membership function tables T2. Referring to this 

ing an internal structure of the memory 12. Referring to figure, the distorted flow model membership function 

the FIG. 4, there are prepared in the memory 12 dis- tables Tl includes K tables T21 to T2K which corre- 
torted flow model membership function tables T2 (see 40 spond to the flow velocity ranges Ui to U*, respec- 

FIG. 7) which are required for allowing the profiles of tively, which are contained in the flow velocity level 

distorted flow distributions to be handled as the fuzzy membership function table Tl, as can be seen in FIG. 4. 

variables, a flow velocity level membership function Shown in FIG. 7 is one of these distorted flow model 

table Tl (FIG. 5) for inferring approximately a flow membership function tables T21 to T2K, e.g. the table 
velocity range of interest, mean flow rate calculating 45 T2k for the flow velocity level Ujt (k is an integer in a 

expression tables T3 (FIG. 10) for calculating the mean range given by 1 <k<K). In the table shown in FIG. 7, 

flow rate in correspondence to the distorted flow model symbols fti, ft* . • • * fttf in the topmost row represent 

pattern determined through the inference (FIG. 10) and normalized values of the flow velocity measurement 

a mean flow rate inference memory part Al for storing values detected by the sensors 51, 52, ... , 5M, respec- 
fuzzy inference procedures. In more concrete, the mean 50 tively. On the other hand, symbols D m t» Dm2. • • * 

flow rate inference memory part Al stores therein the D m * in the leftmost column represent distorted flow 

conditional decision statements and the inference ex- model patterns used for classification of the distorted 

pressions (1) to (7) described hereinbefore. flow distribution profiles. Further, symbols pyk (where 

FIG. 5 shows a structure of the flow velocity mem- i= 1, . . . , N and j= 1, . . . , M) represent the membership 
bership function table Tl. Referring to the figure, sym- 55 function values used for determining the membership 

bois Urfi, Urf2, . . . , VdL entered in the topmost row of grades of M flow velocity measurement values to the 

the table represent, respectively, specified or particular abovementioned distorted flow model patterns, respec- 

flow velocity ranges resulting from classification of the tively. As a concrete example of the method for prepar- 

mean values of the measured flow velocities. On the ing the membership grades ftp, it is conceivable to 
other hand, symbols Ui, U* . . . , Vx in the leftmost 60 describe them in the form of interpolation functions by 

column of the table represent flow velocity levels for making use of the frequency distribution of the flow 

establishing correspondence of the abovementioned velocity measurement values experimentally deter- 

flow velocity ranges to the distorted flow model mem- mined previously or utilize approximate function for 

bership function tables T2. Further, (w(k=l, . . . , K approximating the mean values and variances of the 
and 1= 1, .... L) represent membership function values 65 flow velocity distributions as the representative param- 

for establishing relations between the flow velocity eters by using a Gaussian curve and storing the func- 

levels and the distorted flow model membership func- tions as subroutines. In this case, it is necessary to nor- 

tion tables T2. malize the membership function values by dividing the 
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flow velocity frequency distributions with a maximum 
frequency thereof in conformance with the definition of 
the membership function. 

Now, description will be made of the distorted flow 
model membership function for preparing the distorted 5 
flow model membership function table T2. FIGS. 8A to 
8C show three, typical configurations of the flow path 
which are located upstream of the sensors and which 
provide causes for the generation of plural (three) basic 
distorted flow models. In this case, it is assumed, by way 10 
of example, that three flow sensors 5A, 5B and 5C are 
installed. 

FIGS. 9A to 91 are views for illustrating membership 
functions of the normalized values iu, ujand ucof flow 
velocity measurement values u^, u^and ucoutputted by 15 
the individual sensors 5A, 5B and 5C in response to a 
distorted flow model Dmi generated by a bent pipe 
having a configuration shown in FIG. 8A. Parentheti- 
cally, the distorted flow model is determined in depen- 
dence on the type of fluid of interest, the geometrical 20 
shape of the pipe or conduit, diameter thereof and so 
forth. Of the figures, FIGS. 9A to 9C illustrate member- 
ship functions at a low flow velocity level Ui resulting 
from classification of the flow velocity into, for exam- 
ple, three levels, FIGS. 9D to 9F illustrate membership 25 
functions at an intermediate flow velocity level U2 and 
FIGS. 9G to 91 illustrate membership functions at a 
high flow velocity level U3. As will be seen in these 
figures, when the mean flow rate u corresponds to the 
low flow velocity level U]» the membership value (i.e. 30 
fitness or membership grade) \i\a\ corresponding to the 
normalized value u^oof the flow velocity measurement 
value ua detected by the sensor 5A is determined to be 
\ia from the graph shown in FIG. 9A. In this way, the 
table T2K includes sub-tables containing the member- 35 
ship function values of the normalized flow velocity 
measurement values ui, . . . , um available from the 
outputs of the sensors 51, 52, ... , 5M for the distorted 
flow models Dmi. D m 2, . . . , D m N for a flow velocity 
model Vh The other tables T21 etc. can be prepared in 40 
the similar manner. 

It should here be mentioned that in place of the dis- 
crete membership function values shown in FIGS. 9A 
to 91, linear interpolation functions such as one shown 
in FIG. 9J may be employed. 45 

FIG. 10 shows a structure of the mean flow rate 
calculating expression table T3. This table includes K 
tables T31 to T3k which correspond to the flow veloc- 
ity levels Ui to U* in the table Tl, as can be seen in FIG. 
4. Of these tables, FIG. 10 shows, by way of example, a 50 
table T3k (where k represents an integer in a range 
given by 1 <k<K) for the flow velocity level U*. Each 
of the flow rate calculating expression tables T31 to T3k 
contains the mean flow rate calculating expressions f>jt 
(i= 1, . . . , N and k= 1, • • . , K) for determining the mean 3 ? 
flow rate from the flow velocity measurement values 
(ui, U2, . . . , u„) for the distorted flow model patterns 
Dmi. . . . . D m N at the flow velocity levels Uj, . . . , Ujc, 
respectively. The expression f# may be stored in the 
form of interpolation function by using a polynomial of *® 
the flow velocity measurement values (ui, U2, . . . , mm) 
for the basic distorted flow models at the flow velocity 
levels, respectively, so that the flow rate having the 
maximum distribution frequency can be calculated. 

The expression fa may be given, for example, as fol- 65 
lows: 



where a,fc and b,jt represent coefficients determined for 
every distorted flow model in the tables T31 to T3K. 

Next, referring to FIG. 4, description will be made on 
mutual relations between the flow velocity membership 
function table Tl, the distorted model membership 
function tables T2 and the mean flow rate calculating 
expression tables T3 as well as a data structurizing 
method. The distorted flow model membership func- 
tion table T2 is divided or classified into plural (K) 
function tables T21, T22, . . . , T2K. Similarly, the mean 
flow rate calculating expression table T3 is classified 
into plural (K) mean flow rate calculating expression 

tables T31, T32 T3K. These K function tables T2 

and K expression tables T3 correspond to the flow 
velocity levels Ui, U2, . . . . resulting from division 
of the flow level by K, respectively. In FIG. 4, lines CI, 
C2, . . . , CK represent the abovementioned correspon- 
dence relations, respectively. In the flow velocity level 
membership function table Tl, correspondences are 
established between the measured flow velocity ranges 

Uj, U2, . . . , Ui. and the function tables T21, T22 

T2K and the calculating expression tables T31, T32, . . 
. , T3K through the medium of the weights {*/. Further, 
each of the distorted model membership function tables 
T21, T22, . . . , T2K is divided into plural (N) distorted 
flow model patterns, wherein each of the distorted flow 
model patterns include M membership functions so as to 
be referred to for each of the measurement flow veloc- 
ity values derived from the plural (M) sensor locations. 
Let's consider, for example, the table T2K shown in 
FIG. 4. This table T2K includes plural or N tables 

T2K1, T2K2 T2KN, wherein the table T2KN 

contains plural or M distorted flow model membership 
functions pNmk (where m= 1, . . . , M) with correspon- 
dences being established between each of the mean flow 
rate calculating expression tables T31, T32, . . . , T3K 
and each of the plural (N) distorted flow model pat- 
terns. The table T3K, for example, includes plural or N 
mean flow rate calculating expressions f/* (where i= 1, . 
..,N). 

The mean flow rate inference part Al is designed to 
infer the mean flow rate of a fluid flowing through a 
pipe in accordance with inference procedure incorpo- 
rated therein by referring to the flow velocity level 
membership function table Tl, the distorted flow model 
membership function tables T2 and the mean flow rate 
calculating expression tables T3 structured as described 
above. In FIG. 4, an arrow RT1 indicates that the mean 
flow rate inference part Al refers to the flow velocity 
level membership function table Tl, and arrows Rl, R2, 
. . . , RK indicate that the mean flow rate inference part 
Al refers to the distorted flow model membership func- 
tion tables T2 and the mean flow rate calculating ex- 
pression tables T3. 

FIG. 11 is a functional block diagram of a calculating 
or computing system for illustrating a procedure 
through which the mean flow rate inference part Al 
infers the mean flow velocity as well as interdepen- 
dence relations among variables. In the figure, arrows 
indicate the interdependence relations among the vari- 
ables. A notation "A1~+TT' represents that the mean 
flow rate inference part Al refers to the flow velocity 
level membership function table Tl. Similarly, symbols 
"AWT2" and "A1-+T3" represent reference of the 
tables T2 and T3 by the inference part Al, respectively. 
Since description has already been made on the method 
of determining the mean flow rate, the following de- 
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scription will be directed to the arithmetic operation , U* by repeating the steps F7 to F9 and calculates the 

procedure of the mean flow rate inference part Al with mean value uo from the mean values uq* weighted with 

emphasis being put on the manner in which the mean the membership grades {*in accordance with the infer- 

flow rate inference part Al refers to the flow velocity ence equation (7). The mean value uo thus determined is 

level membership function table Tl, the distorted flow 5 supplied to the processor unit 11. 

model membership function tables T2 and the mean At a step Fll, the arithmetic operation of the infer- 

flow rate calculating expression tables T3. For simplifi- ence part All comes to an end. 

cation of the description, the mean flow rate inference Parenthetically, it should be added that in the first 

part Al is simply termed the inference part Al, while embodiment of the invention described above, the num- 

the flow velocity level membership function table Tl, io berof the distorted flow model patterns may be one (i.e. 

the distorted flow model membership function table T2 fj-i), 

and the mean flow rate calculating expression table T3 Next , a second embodiment of the invention will be 

arc simply called as the tables Tl, T2 and T3, respec- described. 

tively. In FIG. 11, symbols F0 to Fll denote step identi- In conjunction with software internal structure 

e 5?' . . _„ « . , 15 shown in FIG. 4, it should be mentioned that the mean 

Now referring to FIG. 1, the inference part Al starts flow rate ^ ^ mfme6 by only thc distorted 

operation at a step FQ. flow model membership function tobies T2 while omit- 

At a step Fl, the inference part Al fetches the flow ti rcfercncc t0 ^ flow vclocity lcvcl membership 

velocity measurement valuesju, . , uj#) from the function ^ T1 whhout nccd for nodifyillg fa mfcr . 

sensors 51, . . , SM^through the I/O circuit 114 and the 20 ^ expressions of the inference part Al shown in FIG. 

mput ternunals 1111 111M (corresponding to the 4 ^ m$ GmmaM t0 ^ hat L is ^ to ^ K 

step »i in in ^ flow velocity level membership function table Tl 

At a step F2 the inference part Al determines the shown fa pjQ s whcrcin 

mean value u of the flow velocity measurement values 

(U], U2, . . . , mm) in accordance with the expression (1) 25 

(refer to the step S2 shown in FIG. 3). fc/ = 1 (ft = /) <*) 

At a step F3, the inference part Al refers to the table ^ i (* ^ n 

Tl to thereby determine K membership grades f */ cor- 
responding to the measured flow rate range U<ti (where FIG. 12 shows a corresponding structure of the flow 
1 represents one of 1 to L and k= 1, .... K) to which the 30 velocity level membership function table Tl. In prac- 
mean value u belongs (corresponding to the step S3 tice, however, upon execution of the mean flow rate 
shown in FIG. 3). inference by the processor unit 11, it is required to re- 

At a step F4, the inference part Al determines the place the conditional decision statement given by the 
measured flow velocity maximum value \x max (u ]f u 2l . . expression (2) contained in the mean flow rate inference 
. , u M ) (corresponding to the step S4 in FIG. 3). 35 part Al by the following statement: 

At a step F5 ( the inference part Al calculates the if (S € U rf *) THEN | execute k-th one of inference 
normalized flow velocity measurement values uj, . . . , expressions (4)| 

Win accordance with u;«u/ii*« (where j= 1, .... M) whereon the fuzzy inference is performed on the dis- 
(corresponding to the step S5 in FIG. 3). torted flow distribution pattern. At this juncture, Um 

At a step F6, the inference part Al refers to the table 40 represents the k-th measured flow velocity range to 
T2 to determine the membership function values u-</* on which the mean flow rate u in one of the measured flow 
the basis of the values ui, fi 2 , . . . , n M for every distorted velocity ranges U*i to U rf * resulting from the division 
flow model patterns V mU D m 2, . . . D mJV at the flow by K belongs, as described previously. Accordingly, in 
velocity levels Ui, U 2 , . - . , U* (where i = 1, . . . , N, j = 1, the subsequent execution of the inference, there will be 
. . . , M and k=l, . . . , K). 45 selected only the distorted flow model membership 

At a step F7, the inference part Al calculates the function table T2k and the mean flow rate calculating 
membership grades <{w for every distorted flow model expression table T3k. However, since {#=1 (when 
patterns D ml , D m2 , . . . , D„A'at a flow velocity level k=l) and £ w =0 (when Ml), the steps F3 and F10 
U* in accordance with shown in FIG. 11 are rendered unnecessary with the 

SO result that uo* determined in the step F9 represents uo. 
M According to the method described above, inference 

&a » pjjkitij) can be executed at an increased speed because it is sufifi- 

7=1 cient to execute only one of the fuzzy inferences for K 

flow velocity levels, although stability or consistency of 
where t=l, •. . . t N, and j= 1, . . . , M (refer to the step 55 the inferred values for the flow velocity levels may be 
S6 in FIG. 3). degraded more or less. 

At a step F8, the inference part Al refers to the table With the second embodiment of the invention de- 
T3 to determine the mean flow rate calculating expres- scribed above, speeding-up of the inference is intended 
sions f# for every distorted flow model pattern D m i, by using the tables Tl to T3 and the inference expres- 
D m 2> . . * » Dmwat a given flow velocity level U* (where 60 sions of the first embodiment as they are. It will be 
i=l N, and k=l, . . . , K). however be appreciated that the tables and the infer- 

At a step F9, the inference part Al determines the ence expressions stored in the memory 12 may initially 
mean value uo* from the expressions f,* weighted with be so arranged as to speed up the inference execution, as 
the membership grade <f>/* in accordance with the infer- described below in conjunction with a third embpdi- 
ence equation (6) (corresponding to the step S7 in FIG. 65 ment of the invention. 

3). . FIG. 13 is a view showing a structure of the memory 

At a step F10, the inference part Al determines the 12 implemented to speed up the inference in accordance 
mean values uo* for all the flow velocity levels Uw . . . with a third embodiment of the invention, wherein 
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contents of the tables Tl, T2 and T3 are same as those 
shown in FIG. 4. Further, FIG. 14 is a flow chart for 
illustrating the inference procedure according to the 
third embodiment. Referring to the figure, the inference 
part Al calculates a mean value u of the flow velocity 5 
measurement values ui, U2, . . . , uj/to thereby determine 
the measured flow velocity range Vdk to which the 
mean flow velocity value u belongs, as described here- 
inbefore in conjunction with the second embodiment 
(step F2). Subsequently, by referring to only the tables 10 
T2k and T3k specified by the flow velocity range U<«as 
determined, the steps F4 to F9 are executed. The value 
uok determined at the step F9 is outputted as the mean or 
average value uo. 

As will now be understood from the foregoing, the 15 
inference methods according to the embodiments of the 
invention do not reside in determining merely a simple 
mean value of the individual sensor outputs but takes 
into account the characteristics of the distorted flow 
pattern prevailing at the sensor locations upon deter- 20 
mining the mean value of the sensor outputs. By virtue 
of this feature, a plurality of the sensors need not neces- 
sarily be located on a same plane transversal to the axis 
of the pipe or within a predetermined range in the direc- 
tion axially or diameterwise of the pipe. In other words, 25 
the individual sensors 51, 52, ... , 5M may be disposed 
at different positions along the flow of a fluid within the 
pipe with heights h which differ from one to another 
sensor, as is illustrated in FIG. 15. Thus, there can be 
increased the degree of freedom of selecting the posi- 30 
tions internally of the pipe where the sensors are to be 
installed. With the disposition shown in FIG. 15, the 
individual sensors are disposed at different positions 
along the fluid flowing direction, whereby a temporal 
mean value of the flow rate can be determined. 35 

In order to speed up the fuzzy inference, it is effective 
to execute the weighted mean determining arithmetic 
operation (fuzzy integration) at the steps F2, F9 and 
F10 in the inference procedure shown in FIG. 11. FIG. 
16 shows a fourth embodiment of the invention imple- 40 
mented to this end. Referring to the figure, a numerical 
integration processor 13 incorporated in the fuzzy arith- 
metic unit 1 for performing a numerical integration is 
provided and connected to the processor unit 11 via the 
data bus 31. 45 

FIG, 17 is a view showing an air suction control 
system of an internal combustion engine to which an 
embodiment of the flow rate measuring system accord- 
ing to the invention is applied. Referring to the figure, a 
suction air flow denoted by a reference symbol SA is 50 
forced to flow through an air filter 503 contained in an 
air cleaner 511, a suction pipe 100, an air flow meter 
102, a throttle valve 103 and an intake manifold 501 to 
be charged into a combustion chamber or chambers of 
the interna] combustion engine. 55 

Disposed in a main flow path of the air flow meter 
102 are flow velocity sensors 51, 52, . . . , 5M which are 
connected to detection circuits 41, 42, ... , 4M, respec- 
tively, whose outputs are in turn connected through a 
noise eliminating Alter 3 and an analogue-to-digital 60 
(A/D) converter 2 to a fuzzy arithmetic unit 1 which is 
implemented in accordance with one of the first to third 
embodiments of the invention described above and by 
which the mean flow rate of air flowing through the 
suction pipe 100 is calculated by taking into account the 65 
distorted flow components generally denoted by 102. 
The throttle valve 9 disposed at a position downstream 
of the air flow meter 102 is interlocked to an accelerator 



pedal (not shown) for controlling the suction or intake 
air flow. Mounted on the throttle body 103 is an idle 
speed control valve 8 which serves for controlling the 
air flow at the time when the throttle valve 9 is fully 
closed. In FIG. 17, thick solid arrows indicate direc- 
tions of fuel flows. The fuel pumped up from a fuel tank 
505 by a fuel injection pump 506 is injected into within 
the manifold 501 through fuel injectors 507 to be mixed 
with the air flow having passed through the air flow 
meter 103, the resulting air-fuel mixture being charged 
into the engine 500. 

A control unit 510 serves for arithmetically determin- 
ing the amount of fuel injected through the injector and 
the opening degree of the idle speed control valve 8 on 
the basis of the output of the fuzzy arithmetic unit 1, an 
angular signal representing rotation of the throttle valve 
9, an output signal of an oxygen concentration sensor 
508 disposed within an exhaust manifold 511 and an 
output signal of an engine rotation speed sensor 509. In 
accordance with the results of the fuzzy arithmetic 
operation, the injectors 507 and the idle speed control 
valve 8 are correspondingly controlled. 

The air suction pipe system of the internal combus- 
tion engine is often of a much intricate structure because 
of little space available for installation of the piping in 
the engine room, which makes it difficult to provide an 
auxiliary air flow path. Consequently, the air flow 
within the suction pipe contains significant disturbances 
or turbulence, fluidal noise and distorted flow compo- 
nents. By incorporating the flow measurement appara- 
tus according to the invention in this kind of suction 
piping, the mean flow rate can be measured consistently 
with significantly improved accuracy regardless of the 
presence of distorted flows and noise, whereby high- 
precision air-fuel ratio control can be realized. 

As will be apparent from the foregoing, the present 
invention provides advantageous effects mentioned 
below. 

(1) Because of use of fuzzy integration for calculation 
of the intra-pipe mean flow rate, the flow rate value as 
determined is not affected by the output of any particu- 
lar sensor, whereby the mean flow rate value can be 
obtained with high consistency and reliability, being 
protected against influence of local flow disturbances. 

(2) Since the membership grade or fitness grade of a 
distorted flow pattern is determined for each of plural 
classified distorted flow model patterns and since ex- 
tents of contribution of the distorted flow model pat- 
terns to the determination of the mean flow rate by 
using the membership grades (fitness grades) as the 
weights, it is possible to calculate appropriately the 
mean flow rate even for those distorted flow distribu- 
tions which are intermediate the distorted flow model 
patterns. 

(3) The membership functions for the distorted flow 
model patterns can be changed or modified by exchang- 
ing software. As a result, the flow rate measurement can 
flexibly be performed for a variety of distorted flows 
which differ from one another in respect to the distribu- 
tion and the position where the disturbance takes place 
with the sensors being fixed in respect to the shape, 
location of installation, characteristics and other hard- 
ware aspects. 

(4) The flow rate measuring system according to the 
invention can profitably be employed in internal com- 
bustion engines of motor vehicles. In that case, since the 
air flow rate can be measured with a high accuracy even 
in an air piping suction system of an intricate structure, 
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the fuel injection amount can be regulated finely in 
dependence on the measured air flow rate. Thus, an 
air-fuel ratio control assuring a very high accuracy can 
be realized. In general, the geometrical factors such as 
shape of the pipe in the air suction system vary in depen- 5 
dence on the types of the motor vehicles. In this con- 
junction, the flow rate measuring system according to 
the invention can be applied to any types of the vehicles 
without need for any modifications including that of 
data contained in the memory. 10 

The teachings of the invention can be applied to flow 
rate measurement of a variety of fluids regardless of 
whether it is liquid or gas. 

I claim: 

1. A flow rate measuring method, in control appara- 15 
tus for controlling an internal combustion engine, of 
determining an intake flow rate of a fluid flowing 
through a suction pipe of said internal combustion en- 
gine by using a fuzzy inference system which is imple- 
mented on the basis of fuzzy knowledge including fuzzy 
rules each described in the form of an IF part and a 
THEN part, and membership functions defining mean- 
ings of statements described in said IF and THEN parts, 
said control apparatus including speed detecting means 
for detecting rotational speed of said internal combus- 
tion engine, fuel injection means for injecting a fuel into 
the intake air flow within said suction pipe and control 
means for regulating said fuel injection amount of said 
fuel injection means on the basis of weighted flow rate ^ 
means value of said fluid and output signal of said speed 
detecting means for thereby controlling the rotational 
speed of said internal combustion engine said method 
comprising steps of: 

measuring flow velocities of said fluid at a plurality of 35 
locations within said suction pipe; 

preparing a table of membership functions indicating 
grades of fitness of said flow velocities to a dis- 
torted flow model distribution; 

preparing a mean flow rate calculating expression ^ 
giving a flow rate corresponding to said distorted 
flow model distribution; 

determining a membership grade of the measured 
flow velocities to said distorted flow model distri- 
bution by referring to said membership function 45 
table through the medium of said measured flow 
velocities; and 

calculating said weighted flow rate mean value given 
by said mean flow rate calculating expression by 
using as a weight said membership grade to said 50 
distorted flow model distribution. 

2. A flow rate measuring method, in control appara- 
tus for controlling an internal combustion engine, of 
determining an intake flow rate of a fluid flowing 
through a suction pipe of said internal combustion en- 55 
gine by using a fuzzy inference system which is imple- 
mented on the basis of fuzzy knowledge including fuzzy 
rules each described in the form of an IF part and a 
THEN part, and membership functions defining mean- 
ings of statements described in said IF and THEN parts, 60 
said control apparatus including, speed detecting means 
for detecting rotational speed of said internal combus- 
tion engine, fuel injection means for injecting a fuel into 
the intake air flow within said suction pipe, and control 
means for regulating said fuel injection amount of said 65 
fuel injection means on the basis of a weighted means of 
flow rates of said fluid and an outut signal of said speed 
detecting means for thereby controlling the rotational 
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speed of said internal combustion engine, said method 
comprising steps of: 
measuring flow velocities of said fluid at a plurality of 

locations within said pipe; 
preparing a membership function table containing 
membership functions indicating grades of fitness 
of said flow velocities to each of plural distorted 
flow model distributions; 
preparing mean flow rate calculating expressions 
giving flow rates corresponding to said distorted 
flow model distributions, respectively; 
determining membership grades of the measured flow 
velocities to each of said distorted flow model 
distributions by referring to said membership func- 
tion table through the medium of said measured 
flow velocities; and 
calculating said weighted mean of said flow rates 
given by said mean flow rate calculating expres- 
sions by using as weights said membership grades 
to said distorted flow model distributions. 
3. A flow rate measuring method, in control appara- 
tus for controlling an internal combustion engine, of 
determining an intake flow rate of a fluid flowing 
through a suction pipe of said internal combustion en- 
gine by using a fuzzy inference system which is imple- 
mented on the basis of fuzzy knowledge including fuzzy 
rules each described in the form of an IF part and a 
THEN part, and membership functions defining mean- 
ings of statements described in said IF and THEN parts, 
said control apparatus including speed detecting means 
for detecting rotational speed of said internal combus- 
tion engine, fuel injection means for injecting a fuel into 
the intake air flow within said suction pipe, and control 
means for regulating said fuel injection amount of said 
fuel injection means on the basis of weighted mean of 
flow rates of said fluid and output signal of said speed 
detecting means for thereby controlling the rotational 
speed of said internal combustion engine, said method 
comprising steps of: 
measuring flow velocities of said fluid at a plurality of 

locations within said pipe; 
preparing a first membership function table indicating 
grades of fitness of a mean value of said flow veloc- 
ities within said pipe to each of plural flow velocity 
levels; 

preparing second membership function tables giving 
grades of fitness of said flow velocities to a dis- 
torted flow model distribution at each of said flow 
velocity levels; 
preparing mean flow rate calculating expressions 
giving flow rates corresponding to said distorted 
flow model distribution at said flow velocity levels, 
respectively; 

determining a mean value of the flow velocities mea- 
sured at said plurality of locations within said pipe 
to thereby determine on the basis of said mean 
value first membership grades of said mean value to 
said flow levels, respectively, by reference to said 
first membership function table through the me- 
dium of said mean value; 
determining a second membership grade of said flow 
velocity values to said distorted flow model distri- 
bution at each of said flow velocity levels by refer- 
ence to said second membership function tables 
through the medium of said measured flow veloc- 
ity values; and 
calculating said weighted mean of said flow 
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rates given by said mean flow rate calculating expres- 
sions for said flow velocity levels by using as 
weights said first and second membership grades. 

4. A flow rate measuring method, in control appara- 
tus for controlling an internal combustion engine, of 5 
determining an intake flow rate of a fluid flowing 
through a suction pipe of said internal combustion by 
using a fuzzy inference system which is implemented on 
the basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 10 
and membership functions defining meanings of state- 
ments described in said IF and THEN parts, said con- 
trol apparatus including speed detecting means for de- 
tecting rotational speed of said internal combustion 
engine, fuel injection means for injecting a fuel into the 15 
intake air flow within said suction pipe, and control 
means for regulating said fuel injection amount of said 
fuel injection means on the basis of a weighted mean of 
flow rates of said fluid and an output signal of said speed 
detecting means for thereby controlling the rotational 
speed of said internal combustion engine, said method 
comprising steps of: 

measuring flow velocities of said fluid at a plurality of 
locations within said suction pipe; ^ 

preparing a first membership function table contain- 
ing membership functions indicating grades of fit- 
ness of a mean value of said flow velocities within 
said suction pipe to plural flow velocity levels, 
respectively; 30 

preparing second membership function tables con- 
taining membership functions indicating grades of 
fitness of said flow velocities to a plurality of dis- 
torted flow model distributions at said flow veloc- 
ity level, respectively; 35 

preparing mean flow rate calculating expressions 
giving flow rates corresponding to said distorted 
flow model distributions at said flow velocity lev- 
els, respectively; 

determining a mean value of the flow velocities mea- 
sured at said plurality of locations within said suc- 
tion pipe to thereby determine on the basis of said 
mean value first membership grades of said mean 
value to said flow levels, respectively, by reference 
to said first membership function table through the 45 
medium of said mean value; 

determining second membership grades or said flow 
velocity values to said distorted flow model distri- 
butions at said flow velocity levels, respectively, 
by reference to said second membership function 50 
table through the medium of said measured flow 
velocity values; and 

calculating said weighted mean of said flow rates 
given by said mean flow rate calculating expres- 
sions for said flow velocity levels, respectively, by 55 
using as weights said first and second membership 
grades. 

5. A flow rate measuring method according to claim 

4, 

wherein the step of calculating the flow rate within 60 
said pipe includes: 

a step of determining weighted flow rate mean values 
in accordance with said mean flow rate calculating 
expression by using as weights said second mem- 
bership grades as determined; and 65 

a step of deriving a weighted mean of said weighted 
flow rate mean values calculated in accordance 
with said mean flow rate calculating expression for 



said flow velocity levels by using as weights said 
first membership grades determined. 

6. A flow rate measuring system, in control apparatus 
for controlling an internal combustion engine, for deter- 
mining an intake flow rate of a fluid flowing through a 
suction pipe of said internal combustion engine by using 
a fuzzy inference system which is implemented on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and THEN part, and 
membership functions defining meanings of statements 
described in said IF and THEN parts, said control appa- 
ratus including speed detecting means for detecting 
rotational speed of said internal combustion engine, fuel 
injection means for injecting a fuel into the intake air 
flow within said suction pipe, and control means for 
regulating said fuel injection amount of said fuel injec- 
tion means on the basis of an weighted mean of flow 
rates of said fluid and an output signal of said speed 
detecting means for thereby controlling the rotational 
speed of said internal combustion engine, said flow rate 
measuring system comprising: 

a plurality of detecting means for measuring flow 
velocities of said fluid at a plurality of locations 
within said suction pipe; 

storage means containing a table of membership func- 
tions giving grades of fitness of said flow velocities 
to a distorted flow model distribution; 

storage means containing a mean flow rate calculat- 
ing expression giving a flow rate corresponding to 
said distorted flow model distribution; 

means for determining a membership grade of the 
measured flow velocities to said distorted flow 
model distribution by referring to said membership 
function table through the medium of said mea- 
sured flow velocity; and 

means for determining a weighted mean value of said 
flow rates given by said mean flow rate calculating 
expression by suing as a weight said membership 
grade to said distorted flow model distribution. 

7. A flow rate measuring system, in control apparatus 
for controlling an internal combustion engine, for deter- 
mining an intake flow rate of a fluid flowing through a 
suction pipe of said internal combustion engine by using 
a fuzzy inference system which is implemented on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 
and membership functions defining meanings of state- 
ments described in said IF and THEN parts, said con- 
trol apparatus including speed detecting means for de- 
tecting rotational speed of said internal combustion 
engine, fuel injection means for injecting a fuel into the 
intake air flow within said suction pipe, and control 
means for regulating said fuel injection amount of said 
fuel injection means on the basis of an weighted mean of 
flow rates of said fluid and an output signal of said speed 
detecting means for thereby controlling the rotational 
speed of said internal combustion engine, said flow rate 
measuring system comprising: 

a plurality of detecting means for measuring flow 
velocities of said fluid at a plurality of locations 
within said suction pipe; 

storage means for storing a first membership function 
table containing membership functions indicating 
grades of fitness of a mean value of flow velocities 
within said suction pipe to a plurality of distorted 
flow velocity levels, respectively; 

storage means for storing second membership func- 
tion tables containing second membership func- 
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tions indicating grades of fitness of said flow veloc* 
ities to a plurality of distorted flow model distribu- 
tions at said flow velocity levels, respectively; 

means for determining membership grades of said 
flow velocity values to said distorted flow model 5 
distributions, respectively, by reference to said 
second membership function tables through the 
medium of said measured flow velocity values; and 

means for calculating a weighted mean of said flow 
rates given by said mean flow rate calculating ex- 10 
pressions for said flow velocity levels, respectively, 
by using as weights said membership grades. 

8. A flow rate measuring system, in control apparatus 
for controlling an internal combustion engine, for deter- 
mining an intake flow rate of a fluid flowing through a IS 
suction pipe of said internal combustion engine by using 

a fuzzy inference system which is implemented on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 
and membership functions defining meanings of state- 20 
ments described in said IF and THEN parts, said con- 
trol apparatus including, speed detecting means for 
detecting rotational speed of said internal combustion 
engine, fuel injection means for injecting a fuel into the 
intake air flow within said suction pipe, and control 25 
means for regulating said fuel injection amount of said 
fuel injection means on the basis of a weighted mean of 
flow rates of said fluid and an output signal of said speed 
detecting means for thereby controlling the rotational 
speed of said internal combustion engine, said flow rate 30 
measuring system comprising: 
a plurality of detecting means for mesuring flow ve- 
locities of said fluid at a plurality of locations 
within expression by using as weights said second 
membership grades said suction pipe; 35 
storage means containing a first membership function 
table indicating grades of fitness of a mean value of 
said flow velocities within said pipe to a distorted 
flow velocity level; 
storage means containing second membership func- 40 
tion tables indicating grades of fitness of said flow 
velocities to distorted flow model distributions at 
said flow velocity levels, respectively; 
storage means containing a mean flow rate calculat- 
ing expressions giving a flow rate corresponding to 45 
said distorted flow model distribution at each of 
said flow velocity levels; 
means for determining a mean value of the flow ve- 
locities measured at said plurality of locations 
within said suction pipe to thereby determine on 50 
the basis of said mean value a first membership 
grade of said mean value to said flow levels, respec- 
tively, by reference to said first membership funo- 
tion table through the medium of said mean value; 
means for determining a second membership grade of 55 
said measured flow velocity values to said distorted 
flow model distributions for said flow velocity 
levels, respectively, by reference to said second 
membership function table through the medium of 
said measured flow velocity values; and 60 
means for calculating said weighted mean of said 
flow rates given by said mean flow rate calculating 
expressions for said flow velocity level, respec- 
tively, by using as weights said first and second 
membership grades, for thereby determining the 65 
flow rate of the fluid flowing within said pipe. 

9. A flow rate measuring system, in control apparatus 
for controlling an internal combustion engine, for deter- 
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mining an intake flow rate of a fluid flowing through a 
suction pipe of said internal combustion engine by using 
a fuzzy inference system which is implemented on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 
and membership functions defining meanings of state- 
ments described in said IF and THEN parts, said con- 
trol apparatus including, speed detecting means for 
detecting rotational speed of said internal combustion 
engine, fuel injection means for injecting a fuel into the 
intake air flow within said suction pipe, and control 
means for regulating said fuel injection amount of said 
fuel injection means on the basis of a weighted mean of 
flow rates of said fluid and an output signal of said speed 
detecting means for thereby controlling the rotational 
speed of said internal combustion engine, said flow rate 
measuring system comprising: 
a plurality of detecting means for measuring flow 
velocities of said fluid at a plurality of locations 
within said suction pipe; 
storage means containing a first membership function 
table indicating grades of fitness of a mean value of 
said flow velocities within said suction pipe to a 
plurality of distorted flow velocity level, respec- 
tively; 

storage means containing a second membership func- 
tion table indicating grades of fitness of said flow 
velocities to a plurality of distorted flow model 
distributions for said flow velocity levels, respec- 
tively; 

storage means continuing a mean flow rate calculat- 
ing expressions for deriving flow rates correspond- 
ing to said distorted flow model distributions at 
said flow velocity levels, respectively; 
means for determining a mean value of the flow ve- 
locities measured at said plurality of locations 
within said suction pipe to thereby determine on 
the basis of said mean value a first membership 
grades of said mean value to said flow levels, re- 
spectively, by reference to said first membership 
function table through the medium of said mean 
value; 

means for determining a second membership grades 
of said flow velocity values to said distorted flow 
model distributions for said flow velocity levels, 
respectively, by reference to said second member- 
ship function table through the medium of said 
measured flow velocity values; and 
means for calculating said weighted mean of said 
flow rates given by said mean flow rate calculating 
expressions, respectively, by using as weights said 
first and second membership grades, for thereby 
determining the flow rate of the fluid within said 
suction pipe. 

10. A flow rate measuring system according to claim 

> 

wherein said means for calculating the flow rate 

within said pipe includes: 
means for determining weighted flow rate mean val- 
ues in accordance with said mean flow rate calcu- 
lating expression by using as weights said second 
membership grades as determined; and 
means for deriving a weighted mean of said weighted 
flow rate mean values calculated in accordance 
with said mean flow rate calculating expression for 
said flow velocity levels by using as a weight said 
first membership grade as determined. 
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11. A flow rate measuring system according to claim 
9, wherein said plurality of detecting means are dis- 
posed in a same plane extending orthogonally to the 
longitudinal axis of said pipe. 

12. A flow rate measuring system according to claim 
9, wherein said plurality of detecting means are dis- 
posed along a direction in which said fluid flows within 
said pipe. 

13. A flow rate measuring method of determining a 
flow rate of a fluid flowing through an intake pipe of an 
internal combustion engine of an internal combustion 
engine by using a flow rate measuring system including 
a fuzzy inference system which is implemented on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 
and membership functions defining meanings of state- 
ments described in said IF and THEN parts, 

said method comprising steps, performed by said 
flow rate measuring system, of: 



10 



15 



measuring flow velocities of said fluid at a plurality of 20 _ am 



calculating a weighted mean of said flow rates given 
by said mean flow rate calculating expressions by 
using as weights said membership grades to said 
distorted flow model distributions; and 
regulating fuel injection amount of fuel injection 
amount of fuel injection means of said internal 
combustion engine on the basis of said weighted 
mean of said flow rates and a signal indicative of 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 
15. A flow rate measuring method of determining a 
flow rate of a fluid flowing through an intake pipe of an 
internal combustion engine of an internal combustion 
engine by using a fuzzy inference system which is im- 
plemented on the basis of fuzzy knowledge including 
fuzzy rules each described in the form of an IF part and 
a THEN part, and membership functions defining 
meanings of statements described in said IF and THEN 



locations within said pipe; 
preparing a table of membership functions indicating 
grades of fitness of said flow velocities to a dis- 
torted flow model distribution; 
preparing a mean flow rate calculating expression 25 
giving a flow rate corresponding to said distorted 
flow model distribution; 
determining a membership grade of the measured 
flow velocities to said distorted flow model distri- 
bution by referring to said membership function 30 
table through the medium of said measured flow 
velocities; 

calculating a weighted flow rate mean value given by 
said mean flow rate calculating expression by using 
as a weight said membership grade to said distorted 35 
flow model distribution; and 
regulating fuel injection amount of fuel injection 
amount of fuel injection means of said internal 
combustion engine on the basis of said weighted 
flow rate mean value and a signal indicative of 40 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 
14. A flow rate measuring method of determining a 
flow rate of a fluid flowing through an intake pipe of an 45 
internal combustion engine of an internal combustion 
engine by using a flow rate measuring system including 
a fuzzy inference system which is implemented on the 
basis of fuzzy knowledge including fuzzy rules each 
described in the form of an IF part and a THEN part, 50 
and membership functions defining meanings of state- 
ments described in said IF and THEN parts, 
said method comprising steps, performed by said 

flow rate measuring system, of: 
measuring flow velocities of said fluid at a plurality of 55 

locations within said pipe; 
preparing a membership function table containing 
membership functions indicating grades of fitness 
of said flow velocities to each of plural distorted 
flow model distributions; 60 
preparing mean flow rate calculating expressions 
giving flow rates corresponding to said distorted 
flow model distributions, respectively; 
determining membership grades of the measured flow 
velocities to each of said distorted flow model 65 
distributions by referring to said membership func- 
tion table through the medium of said measured 
flow velocities; and 



said method comprising steps, performed by said 

flow rate measuring system, of: 
measuring flow- velocities of said fluid at a plurality of 

locations within said pipe; 
preparing a first membership function table indicating 
grades of fitness of a mean value of said flow veloc- 
ities within said pipe to each of plural flow velocity 
levels; 

preparing second membership function tables giving 
grades of fitness of said flow velocities to a dis- 
torted flow model distribution at each of said flow 
velocity levels; 

preparing mean flow rate calculating expressions 
giving flow rates corresponding to said distorted 
flow model distribution at said flow velocity levels, 
respectively; 

determining a mean value of the flow velocities mea- 
sured at said plurality of locations within said pipe 
to thereby determine on the basis of said mean 
value first membership grades of said mean value to 
said flow levels, respectively, by reference to said 
first membership function table through the me- 
dium of said mean value; 

determining a second membership grade of said flow 
velocity values to said distorted flow model distri- 
bution at each of said flow velocity levels by refer- 
ence to said second membership function tables 
through the medium of said measured flow veloc- 
ity values; 

calculating a weighted mean of said flow rates given 
by said mean flow rate calculating expressions for 
said flow velocity levels by using as weights said 
first and second membership grades; and 
regulating fuel injection amount of fuel injection 
amount of fuel injection means of said internal 
combustion engine on the basis of said weighted 
mean of said flow rates and a signal indicative of 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 
16. A flow rate measuring method of determining a 
flow rate of a fluid flowing through an intake pipe of an 
internal combustion engine by using a flow rate measur- 
ing system including fuzzy inference system which is 
implemented on the basis of fuzzy knowledge including 
fuzzy rules each described in the form of an IF part and 
a THEN part, and membership functions defining 
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meanings of statements described in said IF and THEN 
parts, 

said method comprising steps, performed by said 
flow rate measuring system, of: 

measuring flow velocities of said fluid at a plurality of 
locations within said pipe; 

preparing a first membership function table contain- 
ing membership functions indicating grades of fit- 
ness of a mean value of said flow velocities within 
said pipe to plural flow velocity levels, respec- 
tively; 

preparing second membership function tables con* 
taining membership functions indicating grades of 
fitness of said flow velocities to a plurality of dis 
torted flow model distributions at said flow veloc 
ity level, respectively; 

preparing mean flow rate calculating expressions 
giving flow rates corresponding to said distorted 
flow model distributions at said flow velocity lev- 
els, respectively; 

determining a mean value of the flow velocities mea- 
sured at said plurality of locations within said pipe 
to thereby determine on the basis of said mean 
value first membership grades of said mean value to 
said flow levels, respectively, by reference to said 25 
first membership function table through the me- 
dium of said mean value; 

determining second membership grades or said flow 
velocity values to said distorted flow model distri- 
butions at said flow velocity levels, respectively, 30 
by reference to said second membership function 
table through the medium of said measured flow 
• velocity values; 

calculating a weighted mean of said flow rates given 
by said mean flow rate calculating expressions for 35 
said flow velocity levels, respectively, by using as 
weights said first and second membership grades; 
and 

regulating fuel injection amount of fuel injection 
amount of fuel injection means of said internal 40 
combustion engine on the basis of said weighted 
mean of said flow rates and a signal indicative of 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 

17. A flow rate measuring method according to claim 
16, 

wherein the step of calculating the flow rate within 
said pipe includes: 

a step of determining weighted flow rate mean value 50 
in accordance with said mean flow rate calculating 
expression by using as weights said second mem- 
bership grades as determined; and 

a step of deriving a weighted mean of said weighted 
flow rate mean values calculated in accordance 55 
with said mean flow rate calculating expression for 
said flow velocity levels by using as weights said 
first membership grades determined. 

18. A flow rate measuring system for determining a 
flow rate of a fluid flowing through an intake pipe of an 60 
internal combustion engine by using a flow rate measur- 
ing system including fuzzy inference system which is 
implemented on the basis of fuzzy knowledge including 
fuzzy rules each described in the form of an IF part and 

a THEN part, and membership functions defining 65 
meanings of statements described in said IF and THEN 
parts, 

said flow rate measuring system comprising: 



45 



a plurality of detecting means for measuring flow 
velocities of said fluid at a plurality of locations 
within said pipe; 

storage means containing a table of membership func- 
tions giving grades of fitness of said flow velocities 
to a distorted flow model distribution; 

storage means containing a mean flow rate calculat- 
ing expression giving a flow rate corresponding to 
said distorted flow model distribution; 

means for determining a membership grade of the 
measured flow velocities to said distorted flow 
model distribution by referring to said membership 
function table through the medium of said mea- 
sured flow velocity; 

means for determining a weighted mean value of said 
flow rates given by said mean flow rate calculating 
expression by suing as a weight said membership 
grade to said distorted flow model distribution; and 

control means for regulating fuel injection amount of 
fuel injection means of said internal combustion 
engine on the basis of said weighted mean of said 
flow rates mean value and a signal indicative of 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 

19. A flow rate measuring system for determining a 
flow rate of a fluid flowing through an intake pipe of an 
internal combustion engine by using a flow rate measur- 
ing system including fuzzy inference system which is 
implemented on the basis of fuzzy knowledge including 
fuzzy rules each described in the form of an IF part and 
a THEN part, and membership functions defining 
meanings of statements described in said IF and THEN 
parts, 

said flow rate measuring system comprising: 

a plurality of detecting means for measuring flow 
velocities of said fluid at a plurality of locations 
within said pipe; 

storage means for storing a first membership function 
table containing membership function indicating 
grades of fitness of a mean value of said flow veloc- 
ities within said pipe to a plurality of distorted flow 
velocity levels, respectively; 

storage means for storing second membership func- 
tion tables containing second membership func- 
tions indicating grades of fitness of said flow veloc- 
ities to a plurality of distorted flow model distribu- 
tions at said flow velocity levels, respectively; 

means for determining membership grades of said 
flow velocity values to said distorted flow model 
distributions, respectively, by reference to said 
second membership function tables through the 
medium of said measured flow velocity values; 

means for calculating a weighted mean of said flow 
rates given by said mean flow rate calculating ex- 
pressions for said flow velocity levels, respectively, 
by using as weights said membership grades; and 

control means for regulating fuel injection amount of 
fuel injection means of said internal combustion 
engine on the basis of said weighted means of said 
flow rates mean value and a signal indicative of 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 

20. A flow rate measuring system for determining a 
flow rate of a fluid flowing through an intake pipe of an 
internal combustion engine by using a flow rate measur- 
ing system including fuzzy inference system which is 
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implemented on the basis of fuzzy knowledge including 
fuzzy rules each described in the form of an IF part and 
a THEN part, and membership functions defining 
meanings of statements described in said IF and THEN 
parts, 5 
said flow rate measuring system comprising: 
a plurality of detecting means for measuring flow 
velocities of said fluid at a plurality of locations 
within said pipe; 
storage means containing a lint membership function 10 
table indicating grades of fitness of a mean value of 
said flow velocities within said pipe to a distorted 
flow velocity level; 
storage means containing second membership func- 
tion tables indicating grades of fitness of said flow 15 
velocities to distorted flow model distributions at 
said flow velocity levels, respectively; 
storage means containing a mean flow rate calculat- 
ing expressions giving a flow rate corresponding to 
said distorted flow model distribution at each of 
said flow velocity levels; 
means for determining a mean value of the flow ve- 
locities measured at said plurality of locations 
within said pipe to thereby determine on the basis • 
of said mean value a first membership grade of said 
mean value to said flow levels, respectively, by 
reference to said first membership function table 
through the medium of said mean value; 
means for determining a second membership grade of 30 
said measured flow velocity values to said distorted 
flow model distribution for said flow velocity lev- 
els, respectively, by reference to said second mem- 
bership function table through the medium of said 
measured flow velocity values; 35 
means for calculating a weighted mean of said flow 
rates given by said mean flow rate calculating ex- 
pressions for said flow velocity level, respectively, 
by using as weights said first and second member- 
ship grades, for thereby determining the flow rate 40 
of the fluid flowing within said pipe; and 
control means for regulating fuel injection amount of 
fuel injection means of said internal combustion 
engine on the basis of said weighted mean of said 
flow rates mean value and a signal indicative of 45 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 
21. A flow rate measuring system for determining a 
flow rate of a fluid flowing through an intake pipe of an 50 
internal combustion engine by using a flow rate measur- 
ing system including fuzzy inference system which is 
implemented on the basis of fuzzy knowledge including 
fuzzy rules each described in the form of an IF part and 
a THEN part, and membership functions defining 55 
meanings of statements described in said IF and THEN 
parts, 

said flow rate measuring system comprising: 

a plurality of detecting means for measuring flow 

velocities of said fluid at a plurality of locations 60 

within said pipe; 



storage means containing a first membership function 
table indicating grades of fitness of mean value of 
said flow velocities within said pipe to a plurality of 
distorted flow velocity level, respectively; 

storage means containing a second membership func- 
tion table indicating grades of fitness of said flow 
velocities to a plurality of distorted flow model 
distributions for said flow velocity levels, respec- 
tively; 

storage means continuing mean flow rate calculating 
expressions for deriving flow rates corresponding 
to said distorted flow model distributions at said 
flow velocity levels, respectively; 

means for determining a mean value of the flow ve- 
locities measured at said plurality of locations 
within said pipe to thereby determine on the basis 
of said mean value a first membership grades of 
said mean value to said flow levels, respectively, by 
reference to said first membership function table 
through the medium of said mean value; 

means for determining a second membership grades 
of said flow velocity values to said distorted flow 
model distributions for said flow velocity levels, 
respectively, by reference to said second member- 
ship function table through the medium of said 
measured flow velocity values; 

means for calculating a weighted mean of said flow 
rates given by said mean flow rate calculating ex- 
pressions, respectively, by using as weights said 
first and second membership grades, for thereby 

. determining the flow rate of the fluid within said 
pipe; and 

control means for regulating fuel injection amount of 
fuel injection means of said internal combustion 
engine on the basis of said weighted mean of said 
flow rates mean value and a signal indicative of 
rotational speed of said internal combustion engine, 
thereby controlling the rotational speed of said 
internal combustion engine. 

22. A flow rate measuring system according to claim 
21, 

wherein said means for calculating the flow rate 
within said pipe includes: 

means for determining weighted flow rate mean val- 
ues in accordance with said mean flow rate calcu- 
lating expression by using as weights said second 
membership grades as determined; and 

means for deriving a weighted mean of said weighted 
flow rate mean values calculated in accordance 
with said mean flow rate calculating expression for 
said flow velocity levels by using as a weighted 
said first membership grade as determined. 

23. A flow rate measuring system according to claim 
21, wherein said plurality of detecting means are dis- 
posed in a same plane extending orthogonally to the 
longitudinal axis of said pipe. 

24. A flow rate measuring system according to claim 
21, wherein said plurality of detecting means are dis- 
posed along a direction in which said fluid flows within 
said pipe. 
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